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Abstract In some special zones of the high-dimensional coordinate space of few-body 
systems with identical particles, the operation of an element (or a product of elements) 
of the symmetry groups of the Hamiltonian on a quantum state might be equivalent 
to the operation of another element. Making use of the matrix representations of the 
groups, the equivalence leads to a set of homogeneous linear equations imposing on 
the wave functions. When the matrix of these equations is non-degenerate, the wave 
functions will appear as nodal surfaces in these zones. In this case, these zones are 
prohibited. In this paper, tightly bound 4-boson systems with three types of inter- 
action have been studied analytically and numerically. The existence of the universal 
prohibited zones has been revealed, and their decisive effect on the structures of the 
eigenstates is demonstrated. 
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1 Introduction 

Since various systems may be governed by the same fundamental law, universality 
exists in nature. The early indications of universality in quantum mechanic three- 
body systems was discovered before 1970. 1,2,3 The first strong evidence was given 
by Efimov in 1970 on three weakly bound identical bosons. 05] He found that, when 
the scattering length is sufficiently long, a sequence of loosely bound states, the "Efimov 
states", will appear. The properties of these states are governed by a universal law, 
not depend on the dynamic details of the 3-boson systems. In this paper, another kind 
of universality that exists in tightly bound quantum mechanic few-body systems is 
revealed and confirmed numerically. 

The Hamiltonian of identical particles is invariant under the operations of a set 
of symmetry groups G a , G», ■■■ (including the permutation group). 6,7,8,2] Conse- 
quently, the eigenstates {^(X)} are classified according to the representations of these 
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groups, where X denotes a set of coordinates and i is a serial number. Let be an 
element of G a , gp be that of Gg, and 3 denotes a special zone in the high-dimensional 
coordinate space. When X £ 3, the effects of g a and gp might be equivalent so that 
ga&iiX) — gg]Pi(X). For an example, when 3 is the zone of the squares (SQ), g a 
is a rotation about the normal of the SQ by 2tt/4, and gp is a cyclic permutation of 
particles, then g a and gp are equivalent in 3. Making use of the representations of 
groups, the equivalence leads to a set of homogeneous linear equations 

J2[i%M-D?.M]*i>m = o, (Xes), (l) 

i' 

where Df/^ga) are the matrix elements of the representation. When the matrix of this 
set of equations is non-degenerate, the set S'j) (X) must be zero in 3. In this case, 3 
becomes a prohibited zone (PZ) and the wave function appears as an inherent nodal 
surface (INS). [10] Eq. ^ demonstrates that, for each pair of equivalent operations, 
a constraint will be imposed on the eigenstates. Since the matrixes of representations 
are absolutely irrelevant to dynamics, the constraint is universal disregarding the kind 
of systems (nuclear, atomic, or molecular) and the details of dynamic parameters. It 
implies that the states of different systems but belonging to the same set of represen- 
tations will have exactly the same PZ, and their wave functions will have exactly the 
same INS. On the other hand, some zones are important to binding (say, for 4-body 
systems, the zones of equilateral tetrahedron (ET) and SQ). Whether these zones are 
prohibited is crucial to the binding energy and the geometric character of a state. Fur- 
thermore, the number and the locations of the nodal surfaces in a wave function in 
general determine the strength and mode of oscillation. In particular, the existence of 
the INS implies an inherent mode. Thus, the eigenstates would be seriously affected 
by the universal symmetry constraint. 

The decisive effect of the symmetry constraint on the triply and quadruple excited 
intrashell states of atoms has been revealed previouslv.[ll|[T2, 13 14, 15 Accordingly, 
these states can be naturally classified according to their inherent nodal structures. |16l 
117] For 4-boson systems, a number of predictions on the structures and internal modes 
of oscillation have been made previously. [lSlfl^j However, these predictions have not 
yet been confirmed numerically. The present paper generalizes the work of [18,19] in 
the following aspects: (i) Instead of free 4-boson systems, trapped 4-boson systems 
are considered. Thereby a number of tightly bound states can be obtained which are 
necessary for a systematic analysis, (ii) In addition to theoretical analysis, numerical 
calculations have been performed so that the effect of symmetry constraint can be 
appraised quantitatively, (iii) Three types of interactions have been adopted. The aim 
is to demonstrate the similarity among different kinds of systems. 

In the next section, the symmetry constraints imposing on 4-boson systems are 
studied theoretically. We have chosen appropriate sets of coordinates so that Eq. {T} 
appears in very simple forms and the analysis becomes transparent. Then, an isotropic 
trap together with three types of interaction are introduced, and numerical calculations 
are performed to diagonalize the Hamiltonian. Under the trap the total orbital angular 
momentum L, its Z-component M, and the parity 77 are good quantum numbers. 
Accordingly, an eigenstate with the cm. motion removed can be denoted as iP^M i> 
where i denotes the i-th state of a L^-series. Mostly the i = 1 states (the lowest one) are 
concerned. Therefore, the label i is dropped hereafter (except specified). After obtaining 
the eigenenergies and the eigenstates, a number of quantities (the root-mean-square 
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radius, the one-body densities for the particle distribution, and the shape-densities) are 
further calculated. Thereby, inherent physics can be extracted, and a clear comparison 
among different kinds of system can be made. The emphasis is placed to demonstrate 
the universality of the PZ and the similarity among different systems. A short discussion 
on 4-fermion systems is also given at the end. 

2 Universal prohibited zones in the coordinate space 

The equilateral tetrahedron (ET) and the square (SQ) are the two most important ge- 
ometries. We shall study the symmetry constraint taking place at the ET, SQ, and their 
neighborhoods. The zone associated with the extension-contraction of an ET along one 
of its two-fold (three-fold) axis is defined as Z2-2 (-^1-3)- ^2-2 an d Z1-3 are related to 
the H-type and K-type of oscillations, respectively, that exist in tetramers. [20] When 
the tetrahedrons of Z2-2 ar e further twisted about the two-fold axis, the extended 
zone is defined as Zt w i. 

When appropriate degrees of freedom are chosen, the general expression Eq. |TJ 
could have very simple form. For Z2-2, let the Jacobi coordinates be r = r2 — ri, 
r b = r4 — r 3 , r c = (r 4 + r 3 - r 2 — ri)/2, then r a J- r b _L r c and r a = r b are required. 
r a and r c are considered as variable. We introduce a body frame £' with its origin at 
the cm., its k' lying along r c , j lying along r a , and i' lying along rf,. Z2-2 is shown 
in Fig. QJi. This zone is 2-dimensional in £' . 

In general, an eigenstate can be expressed by using the arguments of a body frame 

as 

*Em(X) = E d% m (WEq(x'), (2) 

Q 

where X denotes the set of coordinates relative to a fixed frame S, and X denotes the 
set relative to £' , Q denotes the set of Euler angles responsible for the transformation 
from £' to [21J Let P denotes an arbitrary particle permutation, 7 a space inversion 
with respect to the cm., and R^ a rotation about the axis v by the angle <f>. In Z2-2 

we have the equivalence IP = ; 2 ■ When IP act at the left side of Eq. (|2| and j 2 
acts at the right side, we obtain the constraint 

H) q *£qPO = n#L Q (x'), (3) 

which is a simple form of Eq. JTJ. Therefore, iP^q^X') is nonzero in Z2-2 om y if 
Q = 0, ±4, • • ■ (when 77 = 1), or Q = ±2, ±6, • ■ • (when i7 = -1). 

We have also the equivalence IP = R^ , which leads to the second constraint 

(-l) L+Q #E t _ Q (X') = IMEq{X'). (4) 

Therefore, the Q — component exists only if (— 1) L = 71. 

With this in mind, we know that ^2-2 is a PZ for 0~ and 1~ states (because they 
have only \Q\ < 1 components), and for 1 + and 3 + states (because they do not have 
\Q\ > 4 components and their Q — component violates the requirement (— 1) L = 77). 

When r c = r a /s/2, the particles form an ET. The zone of ET is denoted as Zet 
which is one-dimensional in E 1 and is a subspace of -£2-2- In ^ET OA (refer to Fig.[l^) 
becomes a three-fold axis of an ET. Thus, additionally, we have 

^2^3 = -^-^27173^ = P ' ( 5 ) 



1 




Fig. 1 (a) The zone Z2-2 defined in a body-frame with r c lying along the third body-axis 
k', r a lying along j', and rj, lying along i'. r a = rj, is assumed. In this zone the magnitudes 
r a and r c are considered as variable, (b) The zone Ztwi defined in a body frame with r c also 
lying along k', r a -L r c lying under the i'-j' plane, and _L r c lying above the i'-j' plane. 
r a = and (f> a = —4>b are assumed. In this zone r c ,r a = ri,, and <j> = (f> a — 4>b are considered 
as variable. 



where /3 = cos 1 ( ^/ 1/3) is the angle between OA and k'. This leads to a set of 
homogeneous linear equations as 

E[(E d Q'Q"(-/ 3 ) e " ¥Q "4"Q(/ 3 ))- 5 Q'Q]^Q'( ET ) = °' ( 6 ) 
Q' Q" 

where <1q,q,, is the matrix elements of the rotation about j'. [21] Under the constraints 
given by Eqs. (|3j), Q and (|6j), nonzero solutions exist at the ET only if L n = + , 3~ 
and 4+ (if L < 4). They are called ET-accessible (ET-ac) states. [X8lfl9] 

When r c — > 0, the shapes in Z2-2 tends to a SQ with its plane lying on the i'-j' 
plane (Fig. |TJi) . The zone of SQ is denoted as -Dsq which is also one-dimensional in S' . 
For -Dsq we have I = P. Thus, the SQ can only exist in II — +1 states. As a subspace 
of the SQ must be constrained by Eqs. © and @. Thus the SQ-ac states are 

0+, 2+ and 4+ (if L < 4). Among them, 0+ and 4+ are both ET- and SQ-ac. 

The geometries in Zi_3 are the regular trihedral pyramids. Let h denotes the height 
which is lying along the 3-fold axis of the pyramid, and s denotes the side length of 
the base regular triangle, h and s are variable. Let k' be lying along the 3-fold axis, 
and j' along a height of the base triangle. Then we have R^ w ^ 3 = P and = IP. 
Accordingly, Zi_3 is a PZ for 0~, 1 + and 2~ states. 

In Z tw ; (refer to Fig. [Tp where r c is lying along k' and can not be seen) r a 1 r c , 
r£, _L r c , r a = rjj are required. r a , r c and the angle <j> between r a and rf, are considered 
as variable. In this zone k' is a 2-fold axis, and i' lying along the bisector of <j> is also 
a 2-fold axis. We have = P and R 1 ^ = P. Accordingly, Zt w j is a PZ for 1~ and 1 + 
states. 

For coplanar shapes, let k' be vertical to the plane of the shapes. Then, R% = I 
must hold. We considered the following three cases. 

(i) The coplanar shape is invariant under inversion (e.g., a parallelogram) and, 
accordingly, k' is a two- fold axis. The associated zone is denoted as Dj±, and the 
equivalence I = P holds. 
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(ii) The shape contains at least a two-fold axis on the i -j plane (say, an isosceles 
trapezoid). The zone is denoted as Dg, and the equivalence = P holds. 

(iii) The intersection of Da and Dg, namely, D^aB- This zone includes the rect- 
angles and diamonds. Specifically, Dsq £ -Dads- 

Associated with the above zones, the PZ are summarized in Table [1] It is em- 
phasized that the existence of the PZ is universal disregarding the kinds of 4-boson 
systems. Since the prohibition is absolutely not violable, the effect of the PZ is decisive. 
The features of all the eigenstates can be thereby more or less predicted. For examples, 
+ is free from symmetry constraint. Thus, it can choose the structure most favorable 
to binding, and its wave function is expected to surround an ET. 2 + is not ET-ac but 
is SQ-ac, therefore its wave function is expected to surround a SQ, and so on. 

For the eight zones listed in Table [T] the most important zones are the Zet an d 
Dgq . The accessibility of these two zones is crucial to the energies of the eigenstates as 
shown below. When a zone is accessible to a state but some subspaces inside the zone 
are not (say, Zt w i is accessible to - but its subspace Z2-2 is not), the zone is not 
a stable zone because the PZ inside the zone implies the existence of inherent nodal 
surfaces, and therefore specific oscillations are excited (examples are given below). 



Table 1 The prohibited zones (PZ) emerge in the coordinate space for the L n states (L < 4) 
of 4-boson systems. The second column gives the dimension of the zones in the body frame. 
The appearance of a PZ is marked by an "/" (say, Z tw ; is a PZ for 1+ and 1 — states). 
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3 Trapped 4-boson systems with three kinds of interaction 

We are going to present numerical results to evaluate quantitatively the consequence 
of the PZ. It is assumed that the four particles are confined by an isotropic harmonic 
trap \rruJ 1 r\. The trap is introduced just for supporting more bound states for a 
systematic analysis, and will not at all affect the appearance of the PZ found above. 
hu) and yjhjmuj are used as units of energy and length. With these units three types of 

interactions are assumed: V A (r) = 10(2e-( r / 1 - 428 ) 2 - e -( r / 2 - 105 ) 2 ), V B {r) = 1000e- 3r - 
40/r 6 (when r > 1.2) and V B (r) = Vb(1.2) (when r < 1.2) , and V c (r) = 15 (when 
r < 1) and Vc(r) = (when r > 1), respectively, where r = |rj — r,-|. Va has a 
short-ranged character and was previously used in nuclear physics for the a-particles, 
Vb belongs to the Van der Waals type for atoms, while Vq is just a repulsive hard core 
potential. In fact, the three interactions are chosen quite arbitrary, just to show the 
possible similarity among different systems. 
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With the interaction and the trap, the Hamiltonian is 

h = E(4 V * + r* 2) + £ Vj{]r * - r ^-> (7) 

i i<j 

where J = A, B, or C. A set of basis functions is introduced to diagonalize the Hamil- 
tonian to obtain the spectra and the eigenstates. The details are given in the Appendix 
A. Note that Va and Vjj both contain a minimum. Thus the total interaction energy 
would be lower if the particle separations are appropriate. Therefore the particles will 
pursue a better geometry. The ET (all six bonds can be optimized) and the SQ (four 
bonds can be optimized) with an appropriate size will be the first and second choices. 
However, for Vq, the minimum is not contained, therefore the pursuit to a better 
geometry is less anxious. 

4 Spectra 

The spectra with the three types of interaction are plotted in Fig. [2] Although the 
interactions are very different, the three spectra have common features. From Table [T] 
we know that the + , 2 + , 4 + , and 3~ are ET-ac and/or SQ-ac. They form the first 
group. This group of states are explicitly lower than the others, and the + state is the 
lowest as expected. Usually, the states with a larger L is higher due to having a stronger 
rotation energy. However, for Va and Vb, 2 + is higher than 3~ because the former is 
simply SQ-ac while the latter is ET-ac. Nonetheless, for Vq, the superiority of ET over 
SQ is weaker, thereby the normal sequence in this group is recovered. For all the three 
spectra, 0~ and 1 + are the highest because many PZ are contained. These two form 
the third group. Although they contain zero or seldom rotation energy, they are still 
much higher than the 4 + . This fact demonstrates the serious effect of the symmetry 
constraint. The remaining four states (1 _ , 2~ , 4~ , and 3 + ) are in the middle, they form 
the second group. Totally speaking, the two spectra for Va and Vb are very similar 
(except the energy scale). The similarity of the spectrum of Vq to the other two is 
weaker. 

5 Root mean square radius 

We calculate the root mean square radius r IBls (L n ) to evaluate the size of each state 
(refer to Appendix B). The ratios r rms are shown in Fig. [3] where the 

three dotted lines guiding the eyes go up and down in a synchronous way. It demon- 
strates once again the similarity among different systems. For all the three types of 
interactions, the sizes of the ET-ac states + , 3~, and 4 + are relatively smaller, while 
the sizes of 0~ and 1 + are the largest (because their wave functions are expelled from 
many PZ as shown below). If dynamic effect played an essential role, the largest would 
have L — 4 because the particles are pushing out by a stronger centrifugal force. But 
in fact not. 

6 Particle distribution 

The particle distributions in the l states are shown by the one-body densities 
r 2 pi(r,9), where r and are the magnitude and the polar angle of a position vector r 
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Fig. 2 Spectra of the 4-boson system with L = to 4 and parity 77 = ±1. Only the lowest 
one of a L^-series are given. Three types of interactions V\, Vb, and Vc are considered. The 
levels have been shifted so that the three ground states (0 + ) of the three types have energy 
zero. The energy unit for Va is !iui. The energy units for the spectra of Vb and Vc have been 
redefined so that the three - states have the same value of energy (but in different units). In 
all the three spectra, the lowest four belong to the first group, the highest two belong to the 
third group. The four in the middle belong to the second group. The states contained in the 
groups are the same for the three types without exception. In a group the same state is given 
by the same kind of line for the three types (say, all the 3~ states belonging to the first group 
are given by the dash-line). 



originating from the cm. (refer to Appendix B). They are shown in Figs. [4] and [5] In 
general the distribution depends on M. The choice M = L implies that the Z-axis is 
chosen to be lying along the direction of L. In Fig. 2] + is isotropic as expected. 1 + 
is nearly isotropic. 3 + and 4 + have a denser distribution when the particle is closer 
to the equator (9 — > n/2). In general, when L is larger, the system will more prefer 
to increase its moment of inertia so as to reduce the rotation energy. This explains 
why 3 + and 4 + prefer to be a little flattened. On the contrary, the distribution of 2 + 



Fig. 3 The ratio r Ima (L n ) /r Tma (0 + ) for the L n states of the three types. Square is for the 
type Va, circle for Vb> an d triangle for Vc- 



prefers to be closer to the Z-axis. This is not due to dynamics but due to symmetry 
constraint. Since 2 + is SQ-ac as shown in Table[T] it will pursue the SQ (This will be 
further confirmed later). When the S' as shown in Fig. Q]is adopted, k' is chosen to 
be vertical to the SQ . However, the SQ can only appear in Q = component due to 
the constraint Eq. (J3]). In this component L is also nearly vertical to k', and therefore 
L and the SQ are nearly coplanar. Thus, in the choice M — L, SQ and the Z-axis are 
nearly coplanar. This explains why the equator is less preferred by 2 + . Furthermore, 
4b is strikingly similar to 4a, 4c is also similar to 4a but in a less extent. Thus, the 
angular distribution of particle is not seriously affected by the details of interaction. In 
the following, we shall see once and once that the similarity between the cases with Va 
and Vb is high, but weaker with Vc- It has been mentioned that, for Vc, the pursuit to 
a better geometry is less anxious. This cause the difference in the extent of similarity. 

High similarity between the cases Va and Vb emerges also in Fig. [5] for the odd- 
parity states. 



7 Shape-densities 

The similarity can be more clearly revealed by observing the wave functions directly. 



Ma — Mb = 1/2 and fj, c = 1 are the reduced masses. Note that R would remain 
invariant under the transformation between two sets of Jacobi coordinates. Therefore 
R is suitable for measuring the size. From the normality and from Eq. ((2j, integrating 
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Fig. 4 One-body densities r 2 pi(r,9) for the even parity states with r being fixed at 

r r ms(L n ). M = L is chosen. Since p(r, 6) = p(r, 7r — 9), only 6 < tt/2 is included. 



over the Euler angles, we have 

1 = (#lm(X)\#lm(X)) = Y, I dX>f Q (X')| 2 . (8) 

There are six arguments included in X' , one is R responsible for the size and the others 
are denoted by S responsible for the shape, and dX' = R s dRdS. Thus the normality 
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Fig. 5 Similar to Fig. [4] but for the odd parity states <]/ 



(-) 



can be rewritten as 1 = J dRdS p(X ), where 



P(X') 



2L + 



T i? 8 5>f Q (X')j 2 , 



(9) 



is the probability density that the system has a given size and a given shape (while 
the orientation of the shape has already been integrated). Therefore p(X') is called the 
shape-density. 
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r 

a 

Fig. 6 The shape-densities p(X') of even-parity states observed in Z2—2 (refer to Fig. |TJi) 
as a function of r a and r c . The three rows are for 0+,2~'~, and 4+. The three columns from 
left to right are for Va_, Vb, and Vq, respectively (the same for Figs.[6]to[9]|. The distribution 
of the density is relative to a body-frame, while the orientation of the body-frame has been 
integrated. The dotted lines have r c = r a /\/2, they denote the ET shapes. In each panel, the 
value of the maximum is given inside a pair of parentheses, while the values of the contours 
form an arithmetic series and decrease to zero. In Fig. 6 - 9 the unit of length is ^/h/mui 
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Fig. 7 The shape-densities p(X') of odd-parity states 2,3, and 4 , respectively, observed 
in Z2—2 as a function of r a and r c . Refer to Fig. [6] 



We shall observe p{X') in some selected zones. Since the ET and the SQ are both 
contained in Z2-2, this zone is important. p(X') of the even-parity and odd-parity 
states observed in Z2-2 are given in Figs. [6] and [3 where r a and r c are variables. Since 
1 + , 3 + , _ , and l - have their wave functions absolutely being zero in this zone, they 
do not appear in these figures (refer to Table[TJ where Z2-2 is a PZ for them). Thus the 
states appear in the figures are + , 2 + , 4 + , 2~ , 3~, and 4~. Note that the + state is 
both ET-ac and SQ-ac. Due to the superiority of the ET over the SQ, all the peaks in 
6a, 6b, and 6c (for Va, Vg, and Vc) are peaked at an ET (where r c = r a /V^) but not a 
SQ (where r c = 0). It has been mentioned that, for Vq, the pursuit to a better geometry 
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is less anxious. Accordingly, the distribution in 6c extends considerably toward a SQ 
. 4 + is also both ET-ac and SQ-ac. The SQ of the 4 + can be nearly vertical to its 
L because |Q| = 4 component is allowed (refer to Eq. ©). Thereby the moment of 
inertia can be increased and the rotation energy can be reduced. Accordingly, for 4 + , 
the SQ becomes more competitive. In fact, for Vc, the SQ has replaced the ET as the 
most probable shape as shown in 6i. 

3~ is ET-ac but not SQ-ac. Accordingly, the distribution is more concentrated 
around the ET and away from the SQ as shown in Figs. [7Ji, [TJ3 , and[7f. 

2 + is not ET-ac but SQ-ac. Accordingly, the distribution is more concentrated 
around the SQ as shown in Figs. [EJl, [EJb, and[5f. 

Among the four states of the second group, 2~ and 4 _ are allowed to appear in 
this zone, and are shown in the first and third rows of Fig. [7] They are neither ET-ac 
nor SQ-ac. There are two peaks in each panel, the peak with a smaller r c is the higher 
peak (except in Fig. [7J) . It implies that the most probable shape is a compressed ET 
along one of its two-fold axis. 

In general, one can define a specific degree of freedom (accordingly, a specific 1- 
dimensional curve in the coordinate space) to describe a specific mode of oscillation. 
When a wave function is remarkably distributed along the curve then the specific mode 
is involved. If a node appears at the curve, the oscillation has been excited. If more 
nodes appear along the curve, the excitation is more vigorous (the simplest example is 
referred to the one-dimensional harmonic oscillator states). In each panel, let a curve 
starting from the upper peak (a prolonged ET), extending to the lower peak, and then 
extending to the abscissa (refer to the solid curve added in Fig. [7g) . Accordingly, r c 
decreases along k (refer to Fig. [TJi) and the prolonged ET suffers a compression and is 
transformed to a SQ. Once r c becomes zero, this process will continue but in reverse 
direction (namely, r c increases further along — k). Then, the SQ is transformed back 
to a prolonged ET, and repeatedly. The mode associated with the transformation is 
named the ET-SQ-ET mode (or the H-type oscillation named in [20]). Since the wave 
functions of 2~ and 4 _ appear as nodes at the ET and SQ, the ET-SQ-ET mode 
contains three nodes and therefore is an excited mode. This explains why 2~ and 4 _ 
are considerably higher than 2 + and 4 + . 

From Table 1 we know that Zi_3 is a PZ for 2~ but not for 4~ . Thus the increase 
of rotation energy in 4~ can be partially compensated due to the Zi_3-ac. Therefore, 
the energy difference between 2~ and 4 _ is small. In fact, the wave function of 4~ 
is also remarkably distributed in ^1-3. Say, for Vc, the maximum of p(X') of 4~ in 
Zi_3 is equal to 1.52 located at a flattened regular pyramid, whereas the maximum in 
Z2-2 is equal to 0.71 located at a prolonged ET (shown in Fig. [7J). Thus, for Vc, the 
most probable shape of 4~ is the flattened pyramid instead of the prolonged ET. 

The other two member 1~ and 3 + of the second group do not appear in ^2-2- 
From the spectra for Va and Vg, we know that 1 _ is higher than 2~ , and 3 + is higher 
than 4~ . This fact implies that the prohibition of Z2-2 will lead to an increase in 
energy. From Table[T]we know that both 1~ and 3 + will be distributed in ^1-3. This 
is shown in Fig. [8^,-|8j; and |8]l-|8f', respectively. In all the panels, there are two peaks. 
The upper peak has larger s (i.e., the base triangle is larger), and is associated with 
a flattened regular pyramid, while the lower peak a prolonged pyramid. Let a curve 
be defined connecting the two peaks. This curve is associated with the transformation 
between the flattened and prolonged pyramid (or the K-type of oscillation [20]). Since 
a node appears at the curve (more exactly, at an ET), this oscillation has been excited. 
Also from Table [1] we know that coplanar structures might exist in these two states. 
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d d d 

1 1 1 

Fig. 8 p{X') of 1 — (a-c) and 3 + (d-f) plotted in Z1—3 for the regular pyramids. The height 
h of the pyramid and the side-length s of the base triangle are variables. The straight lines 
having h = \/2/3s are associated with an ET, and they are an inherent nodal line for these 
states. p(X') of 3+ plotted in Dadb (g-i)i where di and d'2 are the two diagonals of a rhombus. 
The straight line d\ = d2 is associated with the SQ, and it is an inherent nodal line for 3 + . 

This is shown in Figs|5g-|S), where the two peaks are associated with a flattened and a 
prolonged diamond, and a node appears between them. Thus an excited oscillation of 
a diamond around a SQ is contained. 

For the 0~ state of the third group any coplanar structures are prohibited. This 
is due to the equivalence R% = I, where k' is chosen to be vertical to the plane (i?^ 
causes no effect on L = states, while I causes a change of sign on odd parity states, 
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Fig. 9 (a-c) p(X') of _ plotted in a subspace of Z twl (refer to Fig. ^>)> where r c and 
<j> are variables, and r a = depends on r c via a constraint that the hyper-radius R = 

J fia^a + i^b r \ + A*cf c i s fixed at 2r rms (p.s., for an ET, this constraint implies that the distance 

of each particle from the cm. is r rms ). (b) p(X') of 1 + plotted in Da for the parallelograms. 
This zone is a plane. The X- and Y-axes of a body frame with its origin located at the cm. is 
plotted. Two particles are given lying at the X-axis marked by two small circles. The other two 
have coordinates (x, y) and (— x, — y) . The peaks represent the most probable locations of the 
third particle (while the fourth particle at the lower plane can not be seen). Accordingly, the 
most probable shape is a parallelogram. The vertical line at x = represents the diamonds, 
and is an inherent nodal line for 1+. 

this leads to a contradiction) Thus, from Table [T] the wave function of 0~ would be 
essentially distributed in -Ztwi- However, many subspaces in Z tw ; are PZ of 0~. This 
is shown in Figs. [SHUfc. Where <f> = and 7r, and r c — are three nodal lines because 
they are associated with coplanar shapes (refer to Fig.[lp). <j> = 7r/2 is also a nodal line 
because the associated shape belongs also to Z2-2 which is prohibited. Furthermore, let 
the coordinates of the particle 2 relative to £' be denoted as (:r, y, z). Accordingly, those 
of the particles 1, 3, and 4 are (—x,—y,z), (—x,y,~z), and (x,~ y, — z), respectively. 
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One can see that, if y = z the equivalence R^IR 1 ^^ ^ P holds (where P is an 
interchange of a pair of particles), if z = x the equivalence R% IR^, 2 = P holds, and if 

x—y the equivalence R\ IR^ j 2 = P holds. For 0~ , the left side of each of the equivalence 
will cause a change of sign while the right side does not. Therefore 0~ must be zero in 
each of the subspace with y = z, z — x, or x = y. The equality y — z can be rewritten 
as r c /r a = sin(0/2). Note that r c and r a are constrained as y/ r% + r| = 2r rms in 
Figs. HHS=- Thus r c = 2r rm8 sin(</>/2)/V 1 + sin 2 (</>/2) is a nodal line. Similarly, the 
case z=x leads to the nodal line r c = 2r rm s cos (<^/2)/y / l +cos 2 (0/2), while the case 
x=y leads to the nodal line cj> — n/2 mentioned above. Totally, six nodal lines are 
contained in Figs. [9^-[9j;. These lines expel the wave function from favorable zones and 
cause strong inherent oscillation. This explains why the energy of - is so high. 

For 1 + , due to the prohibition of Z tw { and Z1-3, the structure of tetrahedron is 
unfavorable. Furthermore, due to the prohibition of Dg, the coplanar structures are 
also unfavorable. It was found that this state is a mixture of a parallelogram and a 
deformed tetrahedron. The shape-density surrounding the parallelogram is shown in 
Fig.s [5Jl-[5J, where the vertical line is associated with the diamonds and is a nodal line 
for 1 + . The two peaks (each peaked at a parallelogram) separated by the nodal line 
implies an excited oscillation of the parallelogram. 

From the shape-densities plotted in Figs. [6] to [9l the decisive effect of the PZ is 
clearly shown. Since the PZ is universal, the wave functions are expelled in a universal 
way, and the room left for dynamic adjustment is limited. This leads to similarity which 
is clearly shown by comparing the three patterns of any row in the figures for Va, Vg, 
and Vc, respectively. 

8 Fermion systems 

The existence of PZ is not only for 4-boson systems but for all kinds of few-body systems 
containing identical particles. For an example of a 4-fermion system, it is assumed that 
the spin-dependent interaction tends to align the spins (alternatively, a magnetic field 
is applied so that the spins are aligned). For this case, the spatial wave functions will 
be fully anti-symmetric. Accordingly, the equivalent operation Eq. © that hold in 
Z2-2 leads to the fact that $£q(X') is nonzero only if Q = 0, ±4, • ■ • (when II = — 1), 
or Q = ±2, ±6, • ■ • (when 77 = 1). In particular, the Q — component exists only if 
77 = — 1 and L being even. Accordingly, the ET-ac states are now 0~, 3 + , and 4~. 
The SQ-ac states are now 2 + , 3 + , and 4 + , while the states rejected from Z2-2 are 
now + , 1 + , 1~, and 3 _ . Among them 1 + and 3 _ can both access Z1-3, while 1 + 
can access also coplanar structures. Therefore, for the ten states under consideration, 
we can predict that 0~, 2 + , 3 + , and 4~ would be the lowest, while 1~ and + would 
be the highest. We repeat our previous calculation for the case Va but using fully 
anti-symmetric basis functions. The resultant spectra is shown in Table [2] 
The result in Table [2] confirms the prediction. 

9 Final remarks 

From the point of classical dynamics, it is incredible that there are zones in the coordi- 
nate space favorable in energy but absolutely prohibited. However, this paper demon- 
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Table 2 The energies of the ten states (in ascending order) of the fermion system calculated 
with V\ (the unit is hu>). 



State 


o- 


3+ 


2+ 


4- 


1+ 


2" 3- 


4+ 


1- 0+ 


Energy 





0.97 


1.37 


1.77 


1.95 


2.15 2.89 


2.93 


3.22 4.03 



strates analytically and numerically that this is a fact in microscopic few-body systems. 
Since the prohibition arises from the symmetry constraint which depends only on the 
inherent symmetry but is irrelevant to dynamics, the existence of the PZ is universal 
disregarding the kinds of systems and the details of interaction. The quantum states 
of different systems will have exactly the same PZ located at the same subspace when 
these states are governed by the same inherent symmetry (namely, they belong to the 
same set of representations of the symmetry groups). Since some of the PZ have a 
higher geometric symmetry, they are the zones more favorable in energy. Therefore, 
the prohibition causes serious effect and destroys thoroughly the picture painted sim- 
ply by dynamics. Consequently, the structures of the low-lying states are decisively 
determined by their inherent symmetry, and the room left for the adjustment from 
dynamics is narrow, in particular for the i = 1 states 

Although only 4-body systems with [1,1,1,1] and [l 4 ] symmetries are concerned in 
this paper, obviously the existence of the PZ is a universal phenomenon for all kinds 
of few-body systems with identical particles. 

Acknowledgment: The support from the NSFC under the grant No. 10874249 is 
appreciated. 

Appendix 

A Diagonalization of the Hamiltonian 

When the cm. coordinates have been removed, the Hamiltonian for internal motion H\ nt can be 
expressed by a set of Jacobi coordinates r a = r2— ri, r;, = r-4 — r.j, and r c = (r-3+r-4 — ri — r2)/2. 
Let 

MM,r) = -i-V2 + i M r 2 , (10) 
which is the Hamiltonian of harmonic oscillation, Then 

H int = h(l/2,r a ) + h(l/2,r b ) + h(l,r c ) + J2 v (\ r j ( n ) 

i<3 

Let us introduce a variational parameter 7. The eigenstates of ft (7, r) are denoted as <jP, (r) , 
where n and I are the radial and angular quantum numbers, respectively. From (r) a set of 
basis functions for the 4-body system is defined as 

^,^(1234) EE [(^jjra)^^))^^^^)]^, (12) 

where l a and If, are coupled to i (,, then i (, and l c arc coupled to the total orbital angular 
momentum L and its Z-component M, II the parity, k denotes the set of quantum numbers. 



18 



These functions should be further symmetrized. For bosons, l a and ij, should be even and we 
define 

$llILM = X^li7£Af(?>iP2P3P4), (13) 
P 

where the right side is a summation over the permutations. Making use of the Talmi-Moshinsky 
coefficients, 22 , 23 24] each term at the right can be expanded in terms of ^(1234). For an 
example, 

*Z,W1324) = £^*Z',W 1234 ). ( 14 ) 

k' 

where A¥y is the Talmi-Moshinsky coefficients which can be obtained by using the method 
given in the ref. 1241 . Note that the set mM } has not yet been orthonormalized. Thus a 
standard procedure is needed to perform to transform {07 nLM } to a orthonormalized set 

-7 

i&q^n lm} > where q is a serial number to denote a basis function, and each basis function can 

be expanded in terms of {07. ^ iM (1234)}. Finally, the set {0q : nLM} ^ s used to diagonalizc 
-Hint. Due to the Talmi-Moshinsky coefficients, only the set of coordinates r a , r;,, and r c is 
involved. Thus the calculation of the matrix elements is straight forward. When 77 and L are 
given, a series of states will be obtained after the diagonalization. The lowest one of the scries 
is called the first state. 7 has to be adjusted so that the energy of the first state is as low as 
possible. 

Due to having a trap, the bound states will tend to zero quite fast when r tends to infinity. 
For this case, the above basis functions are suitable and will lead to a better convergency. Let 
JV a i,c = 2(n a + rij, + n c ) + la + h + lai and let the number of basis functions be constrained by 
N a b c < N max - Then, as an example of the bosonic case, the resultant energy of the lowest + 
state with the interaction Va is 2.2965, 2.2961, and 2.2960 hw, respectively, when N max = 16, 
18, and 20. Since the emphasis is placed at the qualitative aspect, the convergency appears to 
be satisfactory. 



B One-body density r 2 pi , the root mean square radius r rms 

After the diagonalization, each eigenstate S^jjf is expanded in terms of 0„ nLM> wnere * 
denotes the i-th state of the (77L)-series. For convenience, a transformation from the set 
^q,nLM} to the set {07. ^^^(1234)} is made, and the eigenstate is re-expanded as M = 
Y] k C7,' n ' L 07 nLM (l234). In order to extract the information of a particle from 9^ 'tr, another 
set of Jacobi coordinates r a = r2 — ri , = (ri + T2)/2 — r4, and r e = r*3 — (r4 + ri + r2)/3 
is defined. Accordingly, we have the basis function 

*£W1234) = [W / j (ra)^l 3 ( rd )) !a6 ^l 4 (r e )] iM . (15) 

Since the set {07 ^^^^(1234)} and the set 777^ (1234)} can be transformed into each 
other via the Talmi-Moshinsky coefficients, the eigenstate can once again expanded as 

<ii = J2^ nL Tl nLM (1234), (16) 

k 

where the coefficients DT 1 can be known from the diagonalization and from the transfor- 
mation between sets of basis functions. 
Now, from the normality, we have 

1 = J dr.dr.dr d 0p££)*ip£& 
= J dr sin6»d6» 

[2*(1)V J dr a dr d J2 D2: nL DT nL (T- nLM (123WTl nLM (^)h (17) 
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where 6 is the polar angle of r e , and r = jr e is the distance of particle 3 from the cm. 
Obviously, the quantity inside the bracket is just the one-body density r 2 pi (r, 8) describing 
the particle distribution (all the four particles are distributed in the same way). 
The root mean square radius is 

rrms = 3 { J2 Dll Dll) rW, M M(^)|re 2 |^ iIIM (1234))) 1 / 2 . (18) 
k',k 
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